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Abstract 

The strongly interacting matter under extreme conditions of tem- 
perature for the conjectured deconfined phase of the quark-gluon plasma 
is investigated. A systematic study of the form of the thermal ratio 
of the disorder deviation (TRDD) is presented with the emphasis on 
the degree of deviation of the TRDD-function from unity. The space- 
time evolution of quarks and gluons is studied in the framework of 
two-particle correlation functions to predict the space-time size of the 
closed deconfined phase. 

PACS 12.38.Mh 



1 Introduction 

In the near future a lot of efforts will be needed in searching for the 
quark-gluon plasma (QGP). The existence of a deconfined phase of gluons 
and quarks has been predicted by Quantum Chromodynamics (QCD) [1]. 
Heavy-ion experiments at the CERN Large Hadron Collider (LHC) with the 
program ALICE [2] will start next century, aimed at possible experimen- 
tal confirmation of the hypothetical theoretical predictions of QGP. It is 
supposed that in heavy-ion collisions the hadron matter occurs as a strong 
thermal squeezed state for a short time period. A deconfined phase, QGP, 
is possible at the energy density ~ 2.4 GeV/fm 3 in this thermalized matter 
[3]. However, up to now there are no fixed and guaranteed signals which 
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would make it possible to identify the occurrence of QGP. Matsui and Satz 
[4] pointed out that the signature of the phase transition from the nuclear 
matter to QGP would be the reduction of heavy vector meson yield. There 
is a suppression of a heavy charmonium state J/\l/ in heavy ion collisions via 
the Debye color screening effect in a special mode phase of quarks and gluons- 
the deconfined phase. It is pointed out that charmed quarks and antiquarks 
would leave the zone of deconfinement before forming a hadron composed of 
charmed quarks. 

Following the fruitfull idea of Matsui and Satz [4] we suppose that 

- at multi-TeV (~ 5 — 10 TeV) collider energies the gluon luminosity is far 
larger than the quark luminosity for charmed or beauty hybrid intermediate 
states, where their constituents are distorted temporarily by the deconfined 
environment during the short time scale just after the fusion of two gluons 
coming from colliding ions; 

- the charmed or beauty hybrid mesons turn into the s-wave charmonium or 
bottomonium and 7 quantum. 

Here, in the naive picture, the deconfined state of QGP is occupied by an 
arbitrary number of gluons and/or quarks (g/q). It is clear that there are 
no interactions between gluons and quarks for a short time r if the time 
scale for the QGP existence r = tqgp is larger than the characteristic time 
scale Ty (tqgp > Ty) to form the quark bound state with the average radius 
< r 2 > 1 / 2 =2ry /m q for m q quark mass. Otherwise, there are no unbounded 
quarks if tqgp < 0.32/m and 0.2/m for J/^> and T-particles, respectively. 
The radiative decay of charmed hybrid mesons into J/\l/ with the decay width 
~ 4.0a s keV has been estimated in [5]. There is a very popular point of view 
in literature that a (de) confinement phase transition is predicted to occur at 
the typical energy scale involved, the temperature T c ~ QCD scale Aqcd ~ 
200 MeV. This critical temperature T c is close to the limiting one in hadron 
interactions firstly indicated by Hagedorn [6]. The only remaining problem 
which has been transfering ,Trom one paper to another for a long time is the 
puzzle of rough equality of T c to other three fundamental quantities in QCD: 
the scale Aqcd, the pion mass, and the current mass of a strange quark. 

But, among the issues related to QGP, we attract attention to the problem 
of deconfined phase through the calculation of correlation and distribution 
functions [7] in the thermal theory of quantized fields. In this paper, we 
consider the semiphenomenological model for the QGP existence. To do 
this, we have to use the standard theory of quantized fields replacing: 



1. the asymptotic field operators and 

2. the vacuum expectation values 
by 

1. the thermal field operators and 

2. the thermal statistical averages, 

respectively, in order to formulate correlation and distribution functions of 
produced particles. 

We assume that the heavy-ion collision (ex., Pb-Pb [2]) produces a ther- 
malized quark-gluon gas because of the two-gluon fusion in the local chemical 
and thermodynamic equilibrium. The hot as well as dense gas of particles 
will expand into the surrounding vacuum, thereby cooling down, and the 
particles will interplay with each other and dilute untill free gluons and/or 
quarks, g/ q, as well as particles composed of quarks leave this equilibrium. 
Corrections to the free two-particle correlation functions (CF) taking account 
of the interaction of g/q with their source medium are very actual. This al- 
lows one to get information about the geometry of the emitter-source. For 
convenience, one can distinguish the fase of g/q in the equilibrium from the 
free particles by the freeze-out hypersurface. 

In this paper, we would like to propose new features of formulation with- 
ing the framework of the Langevin-type equation. The method of Langevin 
equation and its extensions to the quantal case have been suggested and 
considered in papers [8-10] and [11-13], respectively. We propose that rather 
complicated real physical processes to happen in the QGP formation should 
be replaced by a one-constituent (ex., gluon or quark) propagation provided 
by a special kernel operator (in the evolution equation) to be considered as 
an input of the model and disturbed by the random force F. We assume 
F to be the external source proposed as both a c-number function and an 
operator. Our evolution equation is an operator one, so that there appear 
new additional problems about the commutation relations and the ordering 
of operators, which do not exist in the classical case [13]. 

Based on the thermal operator-field technique, in Sec. 2 we introduce 
a thermal ratio of the disorder deviation (TRDD), reflecting the degree of 
deviation, from unity, of the ratio of the two-particle thermal momentum- 
dependent distribution to two one-particle thermal distribution functions of 
produced particles, gluons and/or quarks in a partly deconfined phase state. 
We study the four-momentum correlations of two identical particles, g/q, 
which can be both useful and instructive in heavy ion collisions to infer the 



shape of the particle emitter-source (Sec. 3). The sensitivity of the TRDD- 
functions to the size of the emitter will be given in Sec. 3. Within these 
features, the canonical formalism in a stationary state in the thermal equi- 
librium (SSTE) is formulated, and a closed structural resemblance between 
the SSTE and standard quantum field theory is revealed. 



2 Correlation and distribution functions 

To clarify the internal structure of the disordering of particles produced in 
heavy ion collisions, we have to use the consistent approach based on the 
evolution of dynamical variables as well as the extension to different modes 
provided by virtual transitions. 

Let us consider a hypothetical system of the quark-gluon excited local 
thermal phase in QCD where a canonical operator a(k,t) and its Hermitian 
conjugate a + {k,t) occur. In the naive representation , we suppose that in 
this phase there are independent particle sources localized at the points x M 
and all the particles do not interact after their emission. We formulate the 
distribution functions (DF) of produced particles (gluons and quarks) in 
terms of point-to-point equal time temperature-dependent thermal CF of 
two operators 

w(k, k', t- T) = (a + (k, t) a(k', t)) = 

= Tr[a + (k,t) a(k' ,t)e~ H P}/Tr(e- H P) . 

Here, (...) means the procedure of thermal statistical averaging; k and t are, 
respectively, momentum and time variables, e~ Hf3 /Tr(e~ Hl3 ) stands for the 
standard density operator in the equilibrium and the Hamiltonian H is given 
by the squared form of the annihilation a p and creation at operators for 
bose- and fermi- particles, H = J2 P tpO-p^p (the energy e p and operators a p , 
dp carry some index p [14], where p a = 2n n a /L,n a = 0, ±1, ±2, V = L 3 
is the volume of the system considered); (5 is the inverse temperature of the 

— * — * — * — * 

environment, (3 = 1/T. We assume that w(k,k',t;T) — if \k — k'\ < jj, k 
for the characteristic thermolized massive scale where the temperature of 
the environment T does not change. One can suppose that the real physical 
process is devided into two stages: 

1. Thermolization of the hadron matter system or even the quark-gluon 



formation because of the heavy-ion collision. There is an expansion of the 
thermolized space-time as well as a freeze-out process. 

2. Upon cooling to the critical temperature, T c ~ 200°C, the thermolized 
hadron system or the phase with free g/q is decaying into secondary particles 
which could be observed. The first stage, characterizing free particles in 
SSTE, allows one to make any field operator using its expansion like (ex., for 
bose particles) 

where 



tp{xu) = I d k Vk a(k,t) ,vu = =; , (1) 

V " J 1 ' [(2tt) 3 2 A(jfc)]V2 V ; 

and A(k) is the element of the invariant phase volume which will be defined 
later. The standard canonical commutation relation (CCR) 

'a{k,t),a + {k',t)] ± = 5 3 {k-k') (2) 

at every time t is used as usual for bose (-) and fermi (+)-operators. The next 
step is to introduce the TRDD- function D(k, k', t) = R(k, k', t) — 1 reflecting 
the deviation from unity of the R-ratio 

R{k, k', t) = W{k, k', t)/[W(k, t) ■ W{k', t)} . (3) 

The i?-function ( ) is defined as the probability to find two particles, gluons 
or quarks, with momenta k and k! in the same event at the time t normalized 
to the single spectrum of these particles. Here, W(k,t) stands for the one- 
particle thermal DF 

W(k,t) = (b + (k,t)b(k,t)) , (4) 



b(k,t) = a(k,t) + <f)(k,t) (5) 

of the particles emitted (gluons or quarks) under an assumption of occurrence 
of the random source-function <f)(k, t) being an operator, in general. The two- 
particle DF W{k, k',t) looks like 



W(k,k',t) = (b + (k,t) b + (k',t) b(k,t) b(k',t)) 



(6) 



This would allow one to estimate the possibility to find any constituents 
(gluons and quarks) in the excited QCD matter. In order to be applied to 
future experimental search for QGP, it should be important to compare the 
measured TRDD function D for heavy-ion collisions at the TeV scale with 
those calculated in a spherically symmetric model which describes the trans- 
verse momentum spectra of g/ q. Based on the theoretical model considered 
here one can calculate and give the prediction for the space-time size of the 
deconfined system. As a starting point, let us consider, for simplicity, the 
random source-function (p(k, t) in ( g ) as a c-number one. The R-ratio ( |] ) 
of the DF ( H ) and ( f| ) can easily be calculated with the result leading to 

for bose-particles, while for fermi-ones we obtain 
R f (k,k\t) = R b (k,k',t) - 

—* — * 

under the condition a b (k, k',t) > 2cni, where 

c?b(k, k', t) = w(k, k', t) ■ w(k', k, t) + w(k, k', t) ■ <p + (k'; t)<p(k, t) + 
+w(k',k,t) ■ (j) + (k,t)(j)(k',t) , 
a.\ = w(k, k) ■ w(k', k') , «2 = w(k, k) ■ \4>{k')\ , 

« 3 = w(k',k') ■ \m\ 2 ,a 4 = \m\ 2 ■ \<p(k')\ 2 . 

Thus, the deviation from unity, reflecting the naive result, is defined by the 
positive TRDD-functions -Dj,// ranging from to 1: 

n ° h n _ crfe - 2 «i 
u b — —± , u f — • 

Ei=i a* Ej=i «i 

The complete correlation is provided by the condition D b /f=l, while D b /f=0 
corresponds to the naive mode with zeroth correlation. 

The formal structure of the -D&//- functions does not allow one to reveal 
the origin of the particle disordering. We have indicated the formal scheme 
only in order to obtain a nontrivial result (D ^ 0) to reveal the correlation 



effect with the random source represented by the c- number function (p(k;t). 
To avoid the formal representation of the result and to do some translation 
into an operator space, one has to go into some evolution scheme for the set 
of all operators to be used. 

One of our aims is to give a useful tool reflecting the evolution properties 
of propagating particles in a randomly distributed environment. Let us con- 
sider an elementary particle "moving" with the momentum k in the quantum 
equilibrium phase space under the influence of the random force coming from 
surrounding particles. The evolution equations for these particles described 
in terms of the operators b(k, t) and b + (k; t) are 

i d t b(k,t) + A(k,t) = F(k,t) + P . (7) 



- % d t b + (k, t)+A*(k,t) = F + (k,t)+P . (8) 

One can identify both b and b + with the special mode operators of the quark 
and gluon fields or their combinations [15] having the dependence on the 
thermolized QCD matter fields occurring due to heavy-ion collision. In the 
classical picture, equations ( |7| ) and ( || ) are just the Langevin-type ones 
for the Brownian motion of particles. Here, P and F(k,t) stand for the 
stationary external force and the random one, respectively, both acting £from 
the environment. The only operator F has a zeroth value of the statistical 
average, (F) = 0. The correlation between two operators F(k,t) represents 
the fluctuation dissipation theorem of the second kind [16]. The interaction 
of the particles considered with the surrounding ones as well as providing the 
propagation is given by the operator A(k,t) which can be defined as the one 
closely related to the dissipation force. A most simple form of the operator 
A(k,t) should be the following: 

A(k,t)= K(k,t-r) b(k,r) dr . (9) 

J — oo 

Here, an interplay of quarks and gluons with surrounding particles is em- 
bedded into the interaction complex kernel K(k,t), while the real physical 
transitions are provided by the random source operator F(k,t) ( see eq. 
( |7| )). The random evolution field operator K(k,t) in ( [L4| ) stands for the 



random noise and it is assumed to vary stochastically with a 5- like equal 
time correlation function 

(K+(k, r) K(k', t)) = 2 {na) 1/2 k 5{k - k') . (10) 

In the case of the large value of the correlation scale squared a one can get 

(K + (k,r) K(k',r)) -> k exp[-r 2 /(4a)] as a -> oo , (11) 

where = (k — k') while k is the strength of the noise characterized by the 
Gaussian distribution function exp[—^/(4a)]. Hence, both k and a in ( [T0| ) 
and ( |TT] ) define the effect of the Gaussian noise on the evolution of g/q in 
the thermolized environment. 

The formal solution of ( [7] ) in the operator form in <S(3?4) {k^ = (u = 
k°, kj)) is 

~b(k,)=d(k,)+4>(k,) , (12) 

where the operator d(k^) is expressed via the Fourier transformed operator 
F(kft) and the Fourier transformed kernel function K{k lJL ) (coming from ( |9] ) 

as 

a(k lt ) = F(k tt )-[k(k lt )-cj]- 1 , (13) 

while the function ^(k^) is provided by the function ~ P ■ [K{kfj) — u]^ 1 . In 
order to obtain the solution ( [Tj^ ), we have used the fact that at large t the 
terms in the Fourier transformed evolution equation in S(dti), containing the 
fast oscillating factors, weakly tend to zero as t — > +oo as well as t — > — oo. 
The random force operator F(k,t) can be expanded by using the Fourier 
integral 

F(k,t)= / -^(yctye- 1 , 

J-oo Z7T 

where the form ^{k^) ■c(k^) is just the Fourier operator F^k^) = ip(k^) ■c(kfj) 
and the canonical operator c(k^) obeys the commutation relation 

\c(k ll ),c + (k' fl )} ± = 5\k,-k' fl ) . 
The function ^(k^) is determined by the condition 



where A(k) comes from the field expansion ( [I] ), and the CCR ( ) is 

used. The requirement A(k) = (k 2 + m 2 ) in ( |T4[) immediately leads to 
the condition (see also [15]) 



r+°° du 




Loo 2tt 





= 1 . 



The time correlation of the random force operator is defined by 

<f+(m) F(s, t+T )> = r rp d ^-^ 

J-oo J-oo ZTl Z7T 

x(ip + (uj,k) c + (u,k) ip(uj',k) c(u',k)) . 

In the case when the constituents (gluons or quarks) are "moving" inside the 
closed system being in the equilibrium state, the magnitude of the initial time 
t does not have any priority compared to other initial time points. Hence, 
the time correlation p(k,t) = (F + (k,t) F(k,t + r)) should not carry the 
dependence of the time variable. This fact can be realized by putting 

(^ + (lu, k)c + (uj, £)W, k)c(uj', k)) = 2n(c + (k fM )c(k^))\^)\ 2 5(uj - </) . 

Thus, the time correlation is translated into the following one: 

p(k,t) = f +co p- mK)\ 2 <a + (M KK)) ■ 

J —oo Z7T 

In fact, the quantity \ip(k^)\ 2 (c + (k^)c(k fJi )) is the spectral function character- 
izing the random force operator ^{k^). 

3 Distribution functions in S($t±) and space- 
time size. 

In order to study the enhanced probability for emission of two identical parti- 
cles, g/q, and since we are interested in the finite size of the QGP formation, 
the useful information comes from the consideration of DF in S^^). We 
introduce the ratio R of DF as follows: 



(15) 



where W{k li) k'^; T) stands for two-particle probability density subject to 
Bose/Fermi symmetrization and defined as 

W(k„k'^,T) = (b + (k,) ~b + (k'J &(*„) 5(*y> , 

while 

W(k») = (b+(k,) &(*„)) 

is related to the single-particle quantity for a particle with momentum k^. 
Using the Fourier solution of equation ( [7| ) in S(^), one can get R-ratios 
for DF obeying to Bose- 

R b {k^-T) = l + D b {k^-T) (16) 

and Fermi-particles 

^WHW;;T)-2j|^ (it) 



where 

D(k -(k^k'Ji-ik^k,) + 0+^)0^)] + S(A;; ; A;,)0+(^)0(^ 



In order to estimate the thermolized space-time size of QGP we need to have 
the detailed knowledge of CF Using the general approach, one 

cannot do this estimation but suppose the formal expansion to be applied to 
the function E(A;^, k'^), namely 

S(*m,M.) = E<* / IO/(w;T)<C|A;> ) 
c 

where the sum runs over all possible ^-states of the environment while the 
/-function stands for the temperature-dependent one. Using ( [13] ) and the 
Fourier operator = ipik^) ■ c(k^) the two-particle CF 5(/c /i , k') looks 

like 



■ (c + (K) KK)) • (is) 



[k*{k,)-u].[k{k')-u'] v 



The definition of the thermal statistical average (c + (k fl ) c(k')) is related to 
the Kubo-Martin-Schwinger (KMS) condition [17] 

(a{k',t') a + (k,t)) = (a + (k,t) a(k', t - z/3)) • exp(-/3 /i) , (20) 
where fi is the chemical potential. Using the KMS condition ( EDI ) and the 



Fourier transformed solution a{k^) ( [13] ), one can conclude that the thermal 



statistical averages for the c(/c At )-operator should be presented in the following 
form: 

(c + (k fl )c(K)) = 5\k ll -k > tl )-n(u J ,T) , (21) 



(KK) = s\K ~ K) ■ I 1 ± T )] ( 22 ) 

for Bose (+)- and Fermi (-)-statistics where n(u, T) = {exp[(u — //)/3] ±1} -1 . 
Inserting CF ( [L^ ) into ( [L8| ) and taking into account that the 5 4 (k ll — k' )- 
function should be changed by the smooth sharp function Q(r) ■ exp(— g 2 /2), 
one can get the following expression for the -D^-function 

D b {k„, T) = X(k„ T) exp(-g 2 /2) 
x [n{u, T)fi(r) exp(-g 2 /2) + 4>*ihf^%) + 4>*(k,)4>(k'^)] , (23) 

where 

Q(r) 1 
\(ku, k'; T) = -~ K -L n(Cu, T) , Q = -(u + u/) . 

The function O(r) • n(u; T) ■ exp(— q 2 /2) describes the space-time size of the 
QGP fire-ball. Choosing the z-axis along the two-heavy-ion collision axis one 
can put 

q 2 = (r • QoY + (r z ■ Q z f + {r t ■ Q t f , 
Qn = (k — k )^, Qq = tj: — ep, Q z = k z — k z , Qt = \{k x — k x ) + [k y — k y ) ] , 

Q(r) ~ r • r z ■ r\ , 

where t~o, r z and r t are time-like, longitudinal and transverse "size" compo- 
nents of the QGP fire-ball. For estimation the 4-dimensional structure of 
the space-time correlators ( [2T| ), ( |2"2] ) we have to identify r , r z and r t , 



resp., with the averages of the "time" component r, z and \fx 2 + y' 2 on the 
freeze-out hypersurface with the temperature T. Formally, the function 
( ^ ) is the positive one ranging from to 1. Since there are large values of 

\J~Qt = \J (Jhi ~ k'n) 2 m ( El )' a ^ the ^m(^m) dependence should be changed 
by the average value = 0.5(k + k') . The quantitative information (lon- 
gitudinal r z and transverse r t components of the QGP spherical volume, the 
temperature T of the environment) could be extracted by fitting the theoret- 



ical formula ( |23] ) to the measured TRDD function and estimating the errors 
of the fit parameters. Formula ( ^3| ) indicates that a chaotic g/ q source em- 
anating from the thermolized g/q fireball exists. Hence, the measurement of 
the space-time evolution of the g/ q source would provide information of the 
g/q emission process and the general reaction mechanism. In formula ( ^3] ) 
for the Db- function, the temperature of the environment enters through the 
two-particle CF k'^; T). If T is unstable the ^//-functions ( [l^ ) will 

change due to a change of DF W which, in fact, can be considered as an 
effective density of the g/q source. In the case when the temperature in the 
center of the fire-ball is higher than that of the border, this leads to the largest 
contribution to the g/q production. This fact could be "translated into" the 
possible statement that k^ << k z . Formula ( [16] ) looks like the following 
expresion for the experimental R-ratio using a source parametrization: 

R T (r) = 1 + \ T (r) ■ exp(-r t 2 • Q 2 t /2 - r 2 z ■ Q 2 J2) , (24) 

where r t (r z ) is the transverse (longitudinal) radius parameter of the source 
with respect to the beam axis, \t stands for the effective intercept parameter 
(chaoticity parameter) which has a general dependence of the mean momen- 
tum of the observed particle pair. Here, the dependence on the source lifetime 
is omitted. Since < \ T < 1, one can conclude that the effective function At 
can be interpreted as a function of the core particles to all particles produced. 
The chaocity parameter At is the temperature-dependent and the positive 
one defined by 

Mr) n(u;;T)\ 2 

where W(k^) ■ W{k'^) is replaced by Wo^k^, k'^) for convenience regarding the 
point of view that one can distinguish different particles. 



Comparing ( [18] ) and ( |23] ) one can identify 



~ &y = Q{r) ■ n(u>; T) ■ exp(-g 2 /2) . 

In the simple case, the random source function 4>(kfj) in ( [18] ) should be the 
following [18] 



1/2 



(26) 



where a is a real positive parameter. Neglecting <p(k^)- source we immedi- 
ately obtain the trivial result 

D b (k^k'^T) = ~X(uj;T) • exp(-g 2 ) , 

where the temperature-dependent chaoticity X(u;T) is 

n 2 (oo;T) 



A(w;T) 



n(o;; T) • n(u/; T) 



Saving the random source-function we can find the formal representation for 
the TRDD-function Dj, using the a-factorization ( ) 



D b (q 2 ;T) 



A 1/2 (^;T) 
[l + a)(l + a') 



-9 2 /2 



\V 2 (u>;T)e- q2 / 2 + 2(aa') 



1/2 



(27) 



It is easily to see that in the vicinity of q 2 ~ one can get the full correlation 
if a = a' = and A(o>; T)=l. Putting a = a' in ( ) we find the formal lower 
bound on the space-time dimensionless size of the fire-ball for bose-system: 



A>;T) 



[y/(a + l) 2 + « 2 - a] 

In the case of fermi-particles, the following restriction on q 2 is valid (see 
(0)) 



In fact, the function D^k^k'^T) in ( ^3| ) could not be observed in 
the heavy-ion experiment ALICE because of some model uncertainties. In 
the standard consideration, the TRDD-function has to contain a background 
contribution as well as other physical particles (resonances) which have not 
been included in the calculation of the D b - function. In order to be close 
to the experimental data, one has to expand the D b - function as projected 
on some well-defined function (in S^J^)) of the relative momentum of two 
particles produced in heavy-ion collisions 

D^k'^T) - D b (Ql;T) . 

Thus, it will be very instructive to use the polynomial expansion which is 
suitable to avoid any uncertainties as well as characterize the degree of devia- 
tion from the Gaussian distribution, for example. In (—00, +00), a complete 
orthogonal set of functions can be obtained with the help of the Hermite 
polynomials in the Hilbert space of the square integrable functions with the 
measure dfi(z) = exp(—z 2 /2)dz. The function D b corresponds to this class if 

/+00 
dqexp{-q 2 /2) \D b {q)\ n < 00 , n = 0, 1, 2, ... . 
-00 

The expansion in terms of the Hermite polynomials H n (q) 

D b (q) = X J2 Cn • H n (q) ■ exp(-q 2 /2) (28) 

n 

is well suited for the study of possible deviation from both the experimen- 
tal shape and the exact theoretical form of the TRDD function D b ( p3| ). 
The coefficients c n in ( ^8] ) are defined via the integrals over the expanded 
functions D b because of the orthogonality condition 



/+00 
H n {x) H m {x) exp(-x 1 2) dx = 5 n . 
-00 



Hence, the observation of the two-particle correlation (both for bose- and 
fermi-symmetrization) enable to extract the properties of the structure of q 2 , 
i.e. the space-time size of QGP formation. 



4 Conclusion 



1. In this paper we investigated the finite temperature momentum corre- 
lations (of two identical particles, gluons/quarks) which can be both useful 
and instructive in multi-TeV heavy-ion collisions (e.g. Pb+Pb in ALICE 
experimental program) to infer the shape of the gluon/quark source-emitter. 
In fact, we have presented the method of extracting the intercept and source 
parameters from the shape of the TRDD-function. 

2. We used the operator form evolution equation ( |7| ) and ( ^| ) for the special 
mode operators b(k,t) and b + (k,t) (of the gluon/quark field), respectively, 
in the thermolized equilibrium at the freeze-out stage. 

3. The relations between the CF 5(/c^, k' ) and the full i?-functions for bose 
( [16] )- and fermi ( |D] )-particles at the stage of the freeze-out are obtained. 
We have shown the sensitivity of the correlation functions to the space-time 
geometry of the source-emitter ( |2"3| ). In fact, the TRDD- function de- 
scribes the size and shape of the space-time domain where the secondary 
observed particles are generated. 

4. We can conclude that formally, the QGP size scale can be determined 
by the evolution behavior of the field operators and the critical temperature 
T = T c (see formula ( |23| )). 

5. In fact, the full i?-ratio ( [15| ) is the function of four- momentum differ- 
ence of two identical particles as well as the mean total momentum k^. 
Since, the TRDD-function Db is the positive one and restricted by 1, we 
expect that the i?-ratio at too small values of starts from the fixed point 
R{Qli ~~ > 0) = 2 — e (e — > +0) and then falls down (with the Gaussian shape) 
up to unity over some momentum scale interval of an order of the inverse 
source size (see ( |23| ) and ( p5| )). 

6. A large momentum scale determines the minimal size of the QGP 
fire-ball volume for which Rir^T) deviates from its asymptotic value of 
i2(r M ->0;T) = l. 

7. In order to use the quantitative information like longitudinal and trans- 
verse source size, lifetime, we have to fit our expressions to the measured 
R-f unctions. 
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